Current-Induced Membrane Discharge 



M. B. Andersen, 1 ' 2 M. van Soestbergen, 3 * 4 A. Mani, 2 H. Bruus, 1 P. M. Biesheuvel, 3 ' 5 and M. Z. Bazant' 

1 Department of Micro- and Nanotechnology, Technical University of Denmark, 
DTU Nanotech Building 345 East, DK-2800 Kongens Lyngby, Denmark 
2 Department of Mechanical Engineering, Stanford University, Stanford, CA, 94305, USA 
3 Wetsus, Centre of Excellence for Sustainable Water Technology, 
Agora 1, 8934 CJ Leeuwarden, The Netherlands 
4 Department of Applied Physics, Eindhoven University of Technology, 
Den Dolech 2, 5612 AZ Eindhoven, The Netherlands 
5 Department of Environmental Technology, Wageningen University, 
Bornse Weilanden 9, 6708 WG Wageningen, The Netherlands 
6 Departments of Chemical Engineering and Mathematics, 
Massachusetts Institute of Technology, Cambridge, MA, 02139, USA. 
(Dated: March 1, 2012) 

Possible mechanisms for over-limiting current (OLC) through aqueous ion-exchange membranes 
(exceeding diffusion limitation) have been debated for half a century. Flows consistent with electro- 
osmotic instability (EOI) have recently been observed in microfluidic experiments, but the existing 
theory neglects chemical effects and remains to be quantitatively tested. Here, we show that charge 
regulation and water self-ionization can lead to OLC by "current-induced membrane discharge" 
(CIMD), even in the absence of fluid flow. Salt depletion leads to a large electric field which expels 
water co-ions, causing the membrane to discharge and lose its selectivity. Since salt co-ions and water 
ions contribute to OLC, CIMD interferes with electrodialysis (salt counter-ion removal) but could 
be exploited for current-assisted ion exchange and pH control. CIMD also suppresses the extended 
space charge that leads to EOI, so it should be reconsidered in both models and experiments on 
OLC. 

PACS numbers: 47.57.jd, 87.16.dp, 82.45.Mp, 82.33.Ln 



Selective ion transport across charged, water- filled 
membranes plays a major role in ion exchange and desali- 
nation [1, 2], electrophysiology [3], fuel cells [4, 5], and 
lab-on-a-chip devices [6-10], but is not yet fully under- 
stood. A long-standing open question has been to explain 
experimentally observed overlimiting current (OLC), ex- 
ceeding classical diffusion limitation [11]. Possible mech- 
anisms include electroosmotic instability (EOI) and wa- 
ter splitting in the bulk solution [12, 13], as well as sur- 
face conduction and electro-osmotic flow in microchan- 
nels [14] . Vortices consistent with EOI have recently been 
observed under OLC conditions [7, 15, 16], although the 
theory of Rubinstein and Zaltzman [17-19] remains to 
be tested quantitatively. The water splitting mechanism, 
either catalyzed by membrane surface groups or through 
the second Wien effect, has not yet been conclusively tied 
to OLC [13, 20-23]. 

In this Letter, we propose a chemical mechanism for 
OLC, "current-induced membrane discharge" (CIMD), 
resulting from membrane (de)protonation and water self- 
ionization, even in the absence of fluid flow. The ampho- 
teric nature of the charge of ion-exchange membranes 
(i.e. sensitivity to pH and other stimuli) is well known 
[5, 24-29], but not the response to a large applied cur- 
rent. The basic physics of CIMD is illustrated in Fig. 1 
for an anion-exchange membrane. During OLC, a large 
electric field develops on the upstream, salt-depleted side 
of the membrane, which expels H + and attracts OH~, 
causing the membrane to deprotonate and lose selectiv- 
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FIG. 1. [Color online] Basic physics of CIMD, illustrated 
by numerical solutions of Eqs. (2), (5), and (6) for an an- 
ion exchange membrane between two stagnant diffusion lay- 
ers (SDL) for (a) electrostatic potential and concentrations of 
(b) cations c+ and anions c_ and (c) protons ch and hydroxyl 
ions com- 
ity, thereby allowing salt co-ions to pass and producing 
large pH gradients. The upstream solution becomes more 
acidic (low pH), while the downstream, salt-enriched so- 
lution and the membrane become more basic (high pH). 
The local charge of an aqueous membrane strongly de- 
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pcnds on the local pH. In our examples below, we con- 
sider an anion-exchangc membrane with fixed surface 
groups of volumetric concentration c mom , which selec- 
tively allows negatively charged anions (counter-ions) to 
pass, while blocking cations (co-ions) [30]. Depending 
on pH w p[H] = — log 10 (cn), where ch is the proton 
concentration (H + or P±30 + ) in M, the membrane can 
"discharge" (deprotonate): 
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R + H+. 



(1) 



The ratio of product to reactant concentrations in equi- 
librium is the dissociation constant K in M (pK = 
— log 1Q K). Assuming a classical Langmuir adsorption 
isotherm [24-28, 31, 32], the ionization degree of the 
membrane, 



i + — \ = (i + io pH - pK )~ 

chJ 



(2) 



relates its charge concentration a c mcm to pH and pK. 
(For a cation-exchange membrane, the power is pK— pH.) 
To describe the local pH, we cannot assume Boltzmann 
equilibrium with an external reservoir. Instead, we 
consider ion transport coupled to membrane discharge 
Eq. (1) and water self-ionization, 
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with dissociation constant 



K w = c h c h 



(3) 



(4) 



where K w = 1CT 14 M 2 at T = 25 °C. Although kinetics 
can be included [5, 20-22, 33], the reactions (1) and (3) 
are typically fast, so we assume local quasi-equilibrium. 

We now develop a membrane model (seemingly the 
first) including all of these effects: (i) transport of 
four ionic species, including co-ions and water ions (H + 
and OH~) along with majority anions, (ii) water self- 
ionization, and (iii) pH-dcpcndcnt membrane charge. We 
consider the prototypical ID clectrodialysis geometry in 
Fig. 1, consisting of a planar ion-selective membrane of 
thickness L mcm between two well-stirred reservoir com- 
partments of salt ion concentration c rcs and pH of pH ros . 
We adopt the simplest and most commonly used model of 
diffusion limitation [1 1] , in which ion concentrations vary 
across "stagnant diffusion layers" (SDL) of thickness L sc ji 
(of the order 10-100 urn) between the reservoirs and the 
membrane, e.g. representing convection-diffusion bound- 
ary layers or stagnant gel films. 

Ionic diffusion, electromigration and reactions are de- 
scribed by four Nernst-Planck equations. Following 
Rcfs. [24, 25, 34], we combine the Nernst-Planck equa- 
tions for H + and OH~ using Eq. (4) to eliminate the re- 
action terms and relate the water-ion current density J w 
to the water-ion variable c w = (Z?hch — Dqhcoh)/ D w , 



in which D w = \ADh-Doh is the geometric mean of the 
free H + and OH~ diffusivities. We thus arrive at the fol- 
lowing set of coupled, nonlinear, differential equations to 
be solved in both SDLs and the membrane [34]: 
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where Ji is the ionic current density of species i and / r 
is a hindrance factor accounting for porosity, tortuosity 
and constriction (/ r = 1 in the SDLs). Here, <p is the 
dimensionless mean electrostatic potential scaled to the 
thermal voltage Vr = k B T/e = 25.7 mV and satisfying 
Poisson's equation 



d 2 ^ 
dx 2 



-47tA B (Pions + Pmcm) , 



(6) 



where Ab = e 2 /(47re r jEqUbT) is the Bjerrum length, 
and pions = e (c+ - c_ + c H - c h) and p mom = ae c mom 
are charge densities due to the ions and the immobilized 
charges in the membrane, respectively. The porosity e 
of the membrane appears because concentrations Cj are 
defined with respect to the interstitial, not total, volume 
(e = 1 in the SDLs). In our simulations below, we choose 
the following typical parameters: c mcm = 5 M, pK = 9.5, 

imem = ^sdl = 100 UH1, £ riS dl = 78, £ r ,mcm = 29, 

e = 0.4, / r = 0.02 [35], D+ = 1.3 x 10~ 9 m 2 S - 1 
and D_ = 2.0 x 10~ 9 m 2 s _1 (corresponding to NaCl), 
L>h = 9.3 x 10~ 9 m 2 s-\ and D h = 5.3 x 10~ 9 m 2 s~ 1 . 
We also use pH ros = 7 and j3 = 2c ros /c mom = 0.02, un- 
less otherwise noted. The voltage difference across the 
system is A(f>. At the reservoir/SDL boundaries we set 
c± = c rcs and relate c w to pH rcs . 

In spite of neglecting fluid flow, the model still predicts 
OLC, as shown in Fig. 1. The classical ion concentration 
polarization phenomenon is apparent in panel (b) with 
salt depletion where counter-ions (anions) enter (x = X2) 
and enrichment where they leave (x — x 3 ). Within the 
membrane, however, anion depletion and cation (co-ion) 
enrichment reveal a significant loss of selectivity due to 
CIMD. At the same time, panel (c) shows large, order-of- 
magnitudc variations in ch, "mirrored" by coh in equilib- 
rium Eq. (4), with proton enrichment (acidity) in the left 
SDL and proton depletion (basicity) in both the mem- 
brane and the right SDL. The existence of such pH vari- 
ations has been confirmed experimentally in similar sys- 
tems [36-39]. 

Motivated by this observation, we analyze the pH 
gradients perturbativcly in the full CIMD model. We 
consider under-limiting currents, assume thin, quasi- 
equilibrium double layers (Donnan approximation) at the 
SDL/membrane interfaces, and solve the leading-order 
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problem for c+ , c_ and 4> with small perturbations in cjj 
and coh, valid when (ch — coh)/(c+ — c_) <C 1. The 
resulting semi-analytical model (to be described in detail 
elsewhere) suffices to predict CIMD (variations of mem- 
brane charge with local pH) via Eq. (2). Numerical cal- 
culations show that pH and a are nearly constant across 
the membrane, so the water charge density is averaged 
between positions x-i and X3 (see below) to calculate the 
membrane charge and midplane pH [Fig. 2(b)] to be used 
in Eq. (2) to calculate a. 

The final result for the most general model including 
membrane discharge, arbitrary values for pH ros and c ros , 
and the possibility that all diffusion coefficients are dif- 
ferent, consists of Eq. (2) together with the set of al- 
gebraic equations below (see the Supplemental Material 
for details). First, we introduce the dimcnsionless salt 
flux variable j sa it = (J- — J+ D-/D+)/ Ji im , in which 
J\im = — 2 D—C Ies /L s d i is the "classical" limiting current 
density [11], and obtain the salt current- voltage relation, 

A0 = 4tanh- 1 (j salt ) + ^^, (7) 

7 a 

in which 7 = / r // memj where l mcm = L mcm /L sdl 
is the membrane-to-SDL width ratio. The first 
term describes concentration polarization in the SDLs, 
while the second is the Ohmic response of the mem- 
brane. Next, we introduce the dimcnsionless wa- 
ter ion flux j w = J w L s d\ /(fl w V Ky, ) and water 
ion variable p w — c w /\/K w and obtain the follow- 
ing equations, p^xf ™) - p w (x^ cm ) exp[j sa i t (3/(j a)) + 
j w / 7 = 0, sinh- 1 [p w ( 2 ;r m )/2] - sinlr 1 [p w (xf l )/2] - 
sinh- 1 (a/[ ( 0(lTi sa it)]), and p w (xf l ) = p™ T Jw + Po[l + 
27/3/a]ln(l =p Jsait), (where in these expressions i = 2 
and 3 corresponds to — and +, respectively). Here, p™ s is 
related to pH res and p = [4 + (p™ s ) 2 } i . Note that a;™ em 
and x| dl refer to positions on either side of the equilibrium 
electric double layer at the membrane-SDL interfaces. In 
the limit of an infinite membrane charge /3/a — > the 
solution to the leading order problem [Eq. (7)] is simply 
the "classical" result [40], j sa it = tanh(A<^/4). We find 
the characteristic voltage factor <p by expanding Eq. (7) 
for small j sa it <C 1 and obtain j sa i t = A<p/<f>o in which 
4>q = 4 + /3/7 assuming constant a = 1. 

Results of the semi-analytical model are compared 
with full numerical calculations in Fig. 2, which shows 
good agreement in the expected range of validity 
A(f>/4>o < 1. The pH appears to converge towards a 
limiting value for Acf> — > 00, and the jump in this lim- 
iting pH-value between the left SDL and the membrane 
is huge, here about 5 pH units at the highest values of Acj) 
considered. We note that the deviation between the an- 
alytical and numerical solution is largest in the left SDL 
where electroneutrality is most strongly violated. This 
comparative analysis constitutes a validation of our nu- 
merics and provides further support for our conclusions 
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FIG. 2. [Color online] Predicted pH variations from the full 
numerical model, compared to the semi-analytical approxima- 
tion, as a function of the applied voltage (a) in the left SDL, 
just next to the membrane, (b) at the membrane midplane, 
and (c) in the right SDL, next to the membrane. 




FIG. 3. [Color online] Comparison of the classical Mi model 
(only counterions in the membrane) with the full CMID model 
(label "pH" refers to pH rcs ). (a) Membrane ionization degree 
a. (b) co-ion current J+. (c) Total current J to t- (d) Water 
ion current J w . 

regarding the role of pH as controlling the ionic transport 
properties of ion-selective membranes. 

We now turn to a numerical analysis of the CIMD 
model Eqs. (2)-(6). For comparison, we also solve the 
classical model Mi used in all prior work on EOI [15-19] 
in which (i) c_ = in the membrane, (ii) ch = coh = 
everywhere, and (iii) a = 1 for all conditions. We also 
solve two intermediate models which include co-ions in 
the membrane with a = 1 and either exclude (M2) or in- 
clude (M4) water ions, i.e. taking 2 or 4 ions into account 
in the membrane, respectively. The total current density 
is Jtot = J+ + J— + Jw 
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Figure 3(a) shows the significant decrease in the ioniza- 
tion degree a predicted by the CIMD model, in contrast 
to the constant a = 1 in the Mi model. Moreover, a de- 
creases with decreasing /3 (due to increasing Donnan po- 
tential) and decreases with pH ros (due to decreasing ch in 
the membrane). A striking, and yet unexplained predic- 
tion is that for pH ros larger than 7 the ionization degree 
is almost constant until the curve hits that for pH rcs = 7 
after which the curves follow each other. In general we 
find beyond a few times O that reservoir pH has a very 
small influence on membrane charge, fluxes and currents 
(see also Fig. 3(b)-(d)). Figure 3(b) shows the signifi- 
cant increase of co-ion flux J + , thus loss of membrane 
selectivity, with increasing voltage, as predicted by the 
CIMD model, for all values of pH and /3, while Fig. 3(d) 
shows likewise the increase in current density J w due to 
water ions. Still, these contributions do not sum to the 
increased current during OLC, as shown in Fig. 3(c), the 
difference being due to increased counter-ion flux J_ . 

Although the current- voltage relation in CIMD is quite 
complicated, our simulations and analysis suggest two 
general trends: (i) OLC increases with reservoir salt con- 
centration, roughly as /3 65 for the parameters of Fig. 3; 
(ii) OLC is nearly independent of reservoir pH, in spite 
of the large pH gradients produced across the membrane. 

Finally, we analyze the possible effect of CIMD on EOI. 
In the classical Mi model, non-equilibrium space charge 
forms at the limiting current [40-43], and its growing 
separation from the membrane reduces viscous resistance 
to electro-osmotic flow and destabilizes the fluid [17- 
19] . As a measure of the propensity to develop EOI we 
use the transverse (Helmholtz-Smoluchowski) electroos- 
motic mobility /i C o/Moo,o at the left SDL-reservoir edge, 
which is equal to the first moment of the charge density, 
— 47tAb J x 2 x pions dx, or the dimensionless potential dif- 
ference across the left SDL, cj)(xi) — (j>(x2)- 

Figure 4(a) shows that slightly above the limiting cur- 
rent (Jtot/^iim = 1-01) the Mi model already predicts a 
very significant extended space charge layer (the "shoul- 
der" maximum in p iO ns/Pion S ,0 = (Ad L sd i/4> ) d 2 4>/dx 2 
several hundred Debye lengths from the membrane), 
whereas for an even higher current (Jtot/J\im = 1-03), 
using the more realistic CIMD model, the extension of 
this layer is still very minor. The two intermediate mod- 
els lie in between. Figure 4(b) shows how the transverse 
electroosmotic mobility is predicted by the Mi model to 
diverge at the limiting current. This divergence is signif- 
icantly reduced only by the full CIMD model including 
simultaneously co-ion access, water ion transport, wa- 
ter splitting, and membrane discharge. We note that a 
proper analysis of EOI would be more involved, since here 
we have simply focused on the transverse electroosmotic 
mobility as a way of illustrating the suppression of EOI 
due to CIMD. 

In conclusion, we have theoretically demonstrated 
that OLC through aqueous ion-exchange membranes 
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FIG. 4. [Color online] Comparison of three fixed-charge mem- 
brane models M n having n = 1, 2, or 4 mobile ionic species 
and the CIMD model for (a) charge density pions versus dis- 
tance x* from the membrane scaled to the reservoir Debye 
length Ad, and (b) electroosmotic mobility /i C o function of 
total current J to t- 



can result from CIMD, or loss ion selectivity due to 
(de-)protonation coupled to ion transport and water self- 
ionization. The appearance of OLC carried partially 
by salt co-ions and water ions reduces separation effi- 
ciency in electrodialysis, but the associated large pH gra- 
dients and membrane discharge could be exploited for 
current-assisted ion exchange or pH control. CIMD also 
suppresses the non-equilibrium space charge responsible 
for EOI and thus should be considered in both models 
and experiments on OLC with fluid flow. Although we 
have developed the theory for ion-exchange membranes 
in aqueous solutions, CIMD could occur in any nanoflu- 
idic system with an electrolyte whose ions regulate the 
surface charge. 
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